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We numerically investigate an experimentally viable method, that we will refer to as the “chop¬ 
sticks method”, for generating and manipulating on-demand several vortices in a highly oblate 
atomic Bose-Einstein condensate (BEG) in order to initialize complex vortex distributions for stud¬ 
ies of vortex dynamics. The method utilizes moving laser beams (the “chopsticks”) to generate, 
capture and transport vortices inside and outside the BEG. We examine in detail this methodology 
and show a wide parameter range of applicability for the prototypical two-vortex case, and show 
case examples of producing and manipulating several vortices for which there is no net circulation, 
equal numbers of positive and negative circulation vortices, and for which there is one net quantum 
of circulation. We find that the presence of dissipation can help stabilize the pinning of the vortices 
on their respective laser beam pinning sites. Einally, we illustrate how to utilize laser beams as 
repositories that hold large numbers of vortices and how to deposit individual vortices in a sequen¬ 
tial fashion in the repositories in order to construct superfluid flows about the repository beams 
with several quanta of circulation. 


I. INTRODUCTION 

The realm of atomic Bose-Einstein condensates 
(BECs) [IMI has presented a pristine setting where nu¬ 
merous features of the nonlinear dynamics of quantized 
vortices, vortex lattices, and other vortex distributions 
can be both theoretically studied and experimentally ob¬ 
served. Research in this domain has enabled observa¬ 
tions of, for example, precession and excitation of few 
vortices [Mi, collective excitations and dynamics of 
vortex lattices [iiHil, decay of multi ply quantized vor¬ 
tices into singly quantized vortices [TgI, [^ . decay of dark 
solitons into vortices and vortex ring s and gen¬ 

eration of quantum turbulence [2l|, |22|. Additionally, 
experimental efforts directed towards new methods of 
vortex detection HIH are motivated by the need for 
direct measurements of the dynamics of arbitrary dis¬ 
tributions of vortices. These, and numerous other ex¬ 
periments [ 2 ^, demonstrate enormous progress towards 
developing a more complete understanding of vortex dy¬ 
namics in BEGs. 

Nevertheless, to the best of our knowledge, there has 
not been an experimental demonstration of a method to 
construct at will arbitrary (topological charge, as well as 
spatial position) distributions of more than two vortices 
in a BEG. Such a method would enable detailed experi¬ 
mental studies of interactions of vortices with each other, 
with sound, and with trap impurities. New methods to 


*URL: http://nlds.sdsu.edu/ 


study the evolution of many-vortex states involving quan¬ 
tum turbulence, as well as chaotic vortex dynamics [25I - 
[13 would be available, and the role of dissipation in su¬ 
perfluid dynamics could be more precisely determined 
based on experimental data. Here we extend a proof- 
of-principle experimental demonstration, described in a 
companion article [HI, of on-demand vortex generation 
and manipulation of two oppositely charged vortices to 
many-vortex distributions, and further characterize the 
experimental parameters that enable two-vortex gener¬ 
ation and manipulation. We refer to this approach as 
the “chopsticks method”. We consider a highly oblate 
harmonically trapped BEG that is pierced by multiple 
blue-detuned laser beams (that play the role of the chop¬ 
sticks that manipulate the vortices) whose positions and 
intensities can be dynamically controlled. We examine 
conditions for which the motion of the laser beams nu¬ 
cleates vortices in such a way that individual vortices are 
pinned to distinct laser beams during the nucleation pro¬ 
cess. Our numerical results indicate that on-demand en¬ 
gineering of many-vortex distributions is an experimen¬ 
tally realistic possibility, and open up new directions for 
the study of vortex dynamics in BEGs. 

The prototypical case on which our numerical study is 
based involves the presence of two blue-detuned Gaussian 
laser beams that pierce a highly oblate BEG. Each laser 
beam acts as barrier of relatively potential energy Uq that 
of the same order of magnitude than the BEG chemical 
potential —weak beams do not nucleate vortices. Ini¬ 
tially, the two beams are stationary and overlap within 
the BEG. Gonsider a case where both beams begin to 
move in the y direction at a velocity that is significantly 
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lower than the critical speed for vortex dipole nucleation. 
As described in a companion paper , the beams push 
atomic superfluid out of the way, and superfluid fills in 
the space vacated by the laser beam. Simultaneously, 
the beams are pulled apart in the x direction, each beam 
having equal but opposite x-velocity component. Then, 
although the beams are always travelling at a speed well 
below the critical speed for vortex dipole nucleation, the 
“holes” (and the specific path they take) created by the 
laser beams facilitates the formation of two singly quan¬ 
tized vortices of opposite circulation that are simultane¬ 
ously created and pinned (one per beam). Further slow, 
adiabatic, motion (below the critical speed) of the laser 
beams, which serve as vortex optical tweezers, guides the 
position of the pinned vortices. After using the chopstick 
beams to transport the vortices to desired locations, they 
can be released into the BEC to evolve freely by ramp¬ 
ing off the laser beams. It is worthwhile to add here 
that recently an attractive impurity has been utilized in 
a one-dimensional setting, in order to perform similar 
manipulations for one or more dark solitons [29|. 

Here we numerically focus on the use of several laser 
beams to create and manipulate several vortices within 
the BEC. We show that by moving the laser beams to 
desired locations and then ramping them off, releases the 
vortices into the BEC as they subsequently evolve ac¬ 
cording to their inter-vortex dynamics. Since the basic 
vortex generation process creates two vortices of oppo¬ 
site circulation, neutral vortex charge configurations with 
equal numbers of vortices of positive and negative circu¬ 
lation can be readily created. However, driving some of 
the chopsticks beams out of the condensate, non-neutral 
charge configurations can also be generated. Eurther- 
more, by driving multiple chopsticks beams together or 
onto a separate “repository” laser beam, multiply quan¬ 
tized circulations about a single laser beam can be gener¬ 
ated and stored. We confirm that by using either many 
moving laser beams or by depositing vortices to station¬ 
ary repository beams and re-initiating the process of cre¬ 
ating two pinned vortices with two chopsticks beams, ar¬ 
bitrary amounts and configurations of quantized vorticity 
can be prepared in the system, enabling the examination 
of a wide array of associated phenomena. An experi¬ 
mental realization of the prototypical scenario considered 
herein involving two beams is explored in the companion 
paper [^ . 

Our discussion is organized as follows. In Sec. [Ill 
we briefly discuss the theoretical model that is used for 
our study, namely the two-dimensional Gross-Pitaevskii 
equation in the presence of a parabolic trap and a set 
of localized movable defects. In Sec. mil we present our 
numerical results. We start with the simplest case of 
creating two vortices with two chopsticks beams, which 
constitutes our benchmark for quantifying the success of 
the method, and describe how the vortex generation and 
trapping process depends on the beam parameters. Sub¬ 
sequently, we demonstrate that the process can be scaled 
up to neutral configurations of more than two vortices. 


We then consider the removal of a single vortex from neu¬ 
tral configurations by removing one chopstick beam. Eor 
the case of two initial vortices, this allows a single vortex 
to remain in the condensate at a location that is deter¬ 
mined by the remaining beam. More generally, removing 
a single vortex by driving a beam out of the condensate 
leaves an imbalanced vortex charge configuration. Al¬ 
though in the present work, we will consider imbalances 
leading to a total charge of ± 1 , it will be evident that 
the method can enable arbitrary such imbalances in the 
system. Einally, we explore a number of variants to the 
problem. To distill out acoustic energy from these latter 
vortex configurations, we explore the effect of thermally- 
induced dissipation [^, analyzed for vortices, e.g., in 
Ref [sIHsi (see also references therein). We also present 
the possibility of depositing large amounts of vorticity of 
the same charge in a vortex repository. In many cases, we 
explore the vortex dynamics and distributions that result 
from ramping off the laser beams. Einally, in Sec. lEl we 
summarize our findings and discuss directions for future 
study. 

II. MODEL 

We numerically investigate a Bose-Einstein condensate 
in the presence of a strongly anisotropic trapping poten¬ 
tial Wxt = + ^yy‘^ + with ujx = ujy = 

uj±_ uJz- In this case the trapped BEC acquires a nearly 
planar “pancake” shape. We start from the three dimen¬ 
sional Gross-Pitaevskii equation 

■j- 2 

= - —A5' + 53D|4'|"«' + Rext«' (1) 
Zm 

with gsB = , where Og denotes the s-wave scatter¬ 

ing length. The wave function can be factorized into T = 
^( 2 ;)' 0 (x, t), where $( 2 ;) is the ground state of the re¬ 

spective harmonic oscillator: ^(^) = 
and 'ip{x,y,t) is normalized to the number of atoms. Mul¬ 
tiplying Eq. (1) by ^*{z) and integrating over all 2 : yields 
the 2 D GPE: 

indt'ip= - — /\'ip+g2D\i’?i’ + V^^t{x,y,z = Q)'tp. ( 2 ) 

Zm 

Here, we have introduced the 2 D interaction parameter 
^ 2 D = 9zv>jW^Giz) and the harmonic oscillator length 

We non-dimensionalize the 2 D GPE by setting t = tujz^ 
(x^y) = {x^y)/az^ and ^ = Oz^!^. The dimensionless 2D 
GPE then becomes 

+ -v-ip = 0 , (3) 

where g = ^nas/{V^az) and V = ^(x^ where 

the effective ( 2 D) harmonic trapping is Q = u;±/cjz- Fur¬ 
ther letting u = ^/g^^ and dropping all tildes hereafter 
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for notational simplification, the 2 D GPE for the case of 
repulsive interatomic interactions then becomes 

iut + i {Uxx + Uyy) - \u\^ - Vu = Q. (4) 

In the presence of N Gaussian laser beams of identical 
1 /e^ radii a (a dimensionless length, measured in units 
of az)^ the external effective (2D) potential is given by 
the following combination of a harmonic trapping and 
the laser beams: 


V{x,y,t) = 


N 


j = l 


^ x^(t) + y^ (t) 

(5) 


where, for the jth laser beam, Uqj and {xj{t),yj{t)) 
are, respectively, the time-dependent height of the light- 
induced barrier measured in units of hooz and its position 
measured from the trap center in units of az. 

For ease of evaluation of the experimental possibility of 
utilizing our vortex generation and manipulation meth¬ 
ods in the discussion below, we define and utilize a di¬ 
mensional measure of the full width at hal f-maximum of 
the chopsticks beams as = azyj2 ln(2)cr. The beam 
height I/o,i foi* each beam is given in terms of the chem¬ 
ical potential /i measured in units of TiuJz- Furthermore, 
we specify all velocities in terms of the maximum sound 
speed 

C=—V^a^nmax, (6) 

m 

where rimax is the maximum of the three-dimensional 
atomic density (usually attained at the center of the 
parabolic trapping) 

■ ( 7 ) 

Our procedure is conceptually similar to that of the 
companion article [28|, where blue detuned laser beams 
are used in order to produce repulsive barriers that gener¬ 
ate and then dynamically manipulate vortices inside the 
BFG. We specify parameters used in the calculations, 
and present the results of our numerical simulations us¬ 
ing those parameters. 

As mentioned previously, in recent years a significant 
consideration in connection to vortices concerns their dy¬ 
namics in the presence of thermally induced dissipation; 
see, e.g.. Refs, jsol-l^ for which the system is described 
by the dissipative Gross-Pitaevskii equation (dGPE) 

(i - 'j)ut + I (Uxx + Uyy) - \ufu -Vu = 0 , (8) 

where 7 is a dimensionless damping constant. As we 
show below, considering dissipation in our system, not 
only corresponds to a more realistic experimental sce¬ 
nario, but it helps to remove vortices from the periphery 
of the BEG cloud as dissipation induces them to spiral 
outwards. 



FIG. 1: (Color online) Regions of successful and unsuccess¬ 
ful vortex generation and trapping as a function of the beam 
width hyj and the scaled chopsticks speed v/c. The chem¬ 
ical potential corresponds to y = 1.5 and the beam height 
is 2/x. Both beams are initially situated at (0, —20) /xm and 
moved diagonally on straight lines with the corresponding ve¬ 
locity to their final positions at (±10,0) /xm. Region I: The 
velocity of the chopsticks is too small and vortices are not cre¬ 
ated. Region II: Two vortices of opposite charge are created, 
successfully trapped and dragged along with the laser beams. 
Region III: The velocity of the laser beams is too large, result¬ 
ing in the vortices being lost from their respective beams and 
remaining behind the beams. For beam widths approximately 
smaller than bw = ^ /xm no vortex trapping is supported: any 
vortices created are immediately lost from the beams. Movies 
showing the full time evolution exemplary for each of the three 
regimes (and for each of the following figures) can be found 
at http://nonlinear.sdsu.edu/~carreter/Chopsticks.html 


III. NUMERICAL RESULTS 


The initial condition for our evolutionary dynamics 
is obtained by a 2 D fixed point iteration (a Newton’s 
method) in order to identify the ground state of the sys¬ 
tem in the presence of an even number of beams, ranging 
from two to eight. This state is devoid of vortices. Subse¬ 
quently, we compute the time-evolution using a variable- 
order Adams PEGE algorithm, of the type originally 
elaborated in Ref. [sj- 

Although our findings can be straightforwardly gener¬ 
alized to different trapping and atomic gas parameters, 
for concreteness within our pancake-shaped geometry, 
we choose parameter values consonant with the exper¬ 
iments of the companion article [28|. Namely, we choose 
uj± = 27r X 2 Hz, ujz = x 90 Hz, and = 5.3 nm 
corresponding to ^^Rb. We examine a range of chemi¬ 
cal potentials, and indicate this value of /x for each case 
study. These parameters correspond to dimensionless 
times measured in units of 1.77 ms. 
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A. Neutral vortex configurations 

First we discuss the generic example of the creation 
and trapping of a pair of oppositely charged vortices that 
can later be used as a building block for the generation 
of a larger distribution of an even number of vortices, 
with equal numbers of positively and negatively charged 
trapped vortices. Two beams of equal waist and height 
are initially located at the same position and are sub¬ 
sequently moved with the same speed under a suitable 
angle to two different final positions. For a given value 
of the beam size there is a range of suitable values 
of the beam velocity, exceeding a (lower) critical value 
but much less than the speed of sound, for which the 
flow around the beams results in the generation and pin¬ 
ning of two vortices with opposite charge. If the veloc¬ 
ity is chosen to be higher than this range, vortices are 
created but cannot remain trapped by their respective 
beams. These vortices start lagging behind their respec¬ 
tive beams and finally are released from the trapping 
action of the chopsticks beams. The vortices may then 
annihilate one another. Given that this is a prototypi¬ 
cal benchmark scenario, in Fig. [T] we display regions of 
successful and unsuccessful vortex generation and stable 
trapping and manipulation, for various beam widths b^ 
and beam speed v/c [measured in units of the maximum 
speed of sound; see Eq. ([6])] for an experimentally real¬ 
istic chemical potential /i = 1.5, which corresponds to 
^ 4.4 X 10^ atoms comprising the BEG for the parame¬ 
ters given previously, a radial Thomas-Eermi BEG radius 
of ^ 87 yum, and a BEG healing length of ^ 0.64 jam. The 
waist of the beams should be large enough to support the 
existence and trapping of the vortices. We find that for 
the BEG parameters used in our study, with // ^ 1.5, the 
beam width approximately has to exceed 4 fim. Eor beam 
widths below this value, the vortices created are immedi¬ 
ately expelled from the trapping beams [35[ . We have ued 
throughout this study a beam height twice as large as the 
BEG chemical pontental {Uqj = 2/i). However, we have 
checked that for lower beam heights between Uqj = /i 
and Uqj = 2/i, vortex nucleation is still successful for 
beam widths approximately larger than 4 jam. 

We now examine the trajectories of two oppositely 
charged vortices that result from turning off the chop- 
stick beams, in the case of typical symmetric and asym¬ 
metric configurations of the vortices @. We choose a 
value of fi = 1.5 for the chemical potential, a beam 
waist bw = Q jam and a beam height Uq-j = 2ji. It is 
worth mentioning at this stage that lower values of the 
chemical potential resulted in large (density modulation) 
disturbances as the chopsticks move through the conden¬ 
sate, especially if more than two beams are present; on 
the other hand, the case of larger values of the chemi¬ 
cal potential is more difficult to track numerically as the 
vortex width becomes relatively small compared to the 
numerical domain and thus a large numerical grid is nec¬ 
essary. Therefore, for larger numbers of trapped vortices, 
discussed later, it is necessary to increase the value of ji 
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FIG. 2: (Color online) Controlled generation of a vortex 
dipole for ja = 1.5. Throughout this manuscript the upper 
row shows the atomic density in false color (red corresponds 
to highest densities, blue to lowest) while the lower row dis¬ 
plays the phase (blue to red corresponds to phases from 0 to 
27v) at the times indicated in the legend. A pair of beams with 
bw = 6 jam and beam height Uoj = 2/x initially (t = 0) placed 
at (0,-20) jam. are moved within At = 110 to (±10,0) /xm. 
The process nucleates a vortex dipole that is dragged along 
to the desired location. From t = 110 to t = 210 the beams 
are linearly ramped down to release the vortices to the free 
time evolution where the vortices start the typical motion 
of a symmetric dipole pair. The numerics for /x = 1.5 are 
run on 749 x 749 grid points. The resulting spatial discretiza¬ 
tion dx = 0.23 guaranties the existence of about 10 grid points 
within the vortex profile. 



FIG. 3: (Color online) (a) Vortex trajectories for the symmet¬ 
ric vortex configuration from Fig. [2] displayed from the release 
of vortices at t = 210 up to t = 4000. (b) Vortex trajecto¬ 
ries for the asymmetric configuration of Fig. [J] from release of 
vortices up to t = 16000. 


at the expense of more intensive numerics. The typical 
generation of a symmetric vortex configuration is shown 
in Fig. [21 After the vortices have been created and have 
reached the desired symmetric final positions (±10, 0) jam 
the beams are adiabatically and linearly ramped down to 
release the vortices to undergo free (i.e., uninhibited by 
the presence of the chopsticks beams) time evolution. As 
a result, the released vortices exhibit the typical dynam¬ 
ical features of a vortex dipole configuration [ol, [36| . This 
is illustrated in Fig. |3l[a) for the symmetric case of Fig. |2l 
Figured] shows a similar example as the one depicted in 
Fig. [2] but for an asymmetric motion of the chopsticks. 
In this case, one laser beam is kept fixed after t = 110 at 
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FIG. 4: (Color online) Similar to Fig. [2] but for an asymmet¬ 
ric vortex configuration. While one beam is kept fixed after 
110 at (10, 0) /xm, the other beam is moved further until 
t = 200. Then both beams are kept fixed in this asymmetric 
configuration and are linearly ramped down within At = 100. 
As a result, asymmetric vortex dipole dynamics arises past 
the ramp-down time. 


(10,0) /xm while the other beam is moved further until 
t = 200 where both beams are kept fixed and ramped 
down. This procedure seeds an asymmetric configura¬ 
tion that, after removal of the chopsticks, evolves in the 
typical ^itrochoidal trajectories for asymmetric vortex 
dipoles [sl, [sil as shown in Fig[3l^b). 

Next, we investigate the possibilities to create larger, 
even, numbers of vortices in a neutral configurations with 
four, six, and eight vortices. The idea is to start with two, 
three or four pairs of overlapping chopstick beams close 
to the center of the BEG and use the same methodology 
as above to create vortices from each beam pair with the 
same protocol as described above. Doing so, each chop- 
stick nucleates and moves an independent vortex that 
might be placed in a desired location. The beams then 
are kept fixed at their final destination and subsequently 
linearly and adiabatically ramped down to release the 
vortices that in turn are free to evolve without the chop¬ 
sticks being present. Figure [5] shows that this method¬ 
ology is indeed feasible for controllably creating, moving 
and releasing configurations bearing four, six, and eight 
vortices. In principle, this method can be straightfor¬ 
wardly generalized to even larger number of vortices as 
long as there is enough room within the BEG to move the 
chopsticks beams. As it can be observed from the last 
density distributions, the motion and removal (ramping 
down) of multiple chopsticks has a significant perturba¬ 
tive, effect on the background density. In fact, if the 
disturbances to the condensate are considered to be too 
strong, it would be advantageous to increase the value 
of the chemical potential in order to diminish the role of 
interference of the resulting sound waves. Eor this pur¬ 
pose, in our examples with multiple chopsticks we choose 
a slightly larger value of the chemical potential (when 
compared to the previous results with a single chopstick 
pair) of /X = 2. In view of the disturbances created by a 
larger number of beams, we found that starting close to 



FIG. 5: (Golor online) Gontrolled generation of four, six, and 
eight vortices. Ghosen parameters are: chemical potential /x = 
2, beam heights 2/x, and beam widths bw = 6 Mm. Upper set 
of panels: Two pairs of beams initially located at (0, ±8) /j^m 
are moved outward with v/c = 0.123 until t = 80, then their 
positions are kept fixed. From t = 80 until t = 280 the beams 
are linearly ramped down, the vortices are finally released 
and the time-evolution is monitored. Middle panel: Similar 
but for 3 pairs of beams initially located symmetrically at a 
distance of 10 fim from the origin with linear ramp down from 
t = 70 until t = 270. Lower panel: Similar but for 4 pairs of 
beams initially symmetrically located at a distance of 15 fim 
from the origin and moved at a speed x;/c = 0.15 and linearly 
ramped down from t = 50 until t = 250. 


the center of the BEG and moving the beams outward is 
advantageous in comparison to starting further outwards 
and moving the beams inward. 


B. Non-neutral configurations 

We now consider creating an odd number of vortices 
with a charge imbalance of one, such that there is one 
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FIG. 6 : (Color online) Controlled generation of a single vortex 
for /i = 2, beam height 2/i and beam width bw = 6 
The initial pair of beams located at (0, —20) fim is moved 
with v/c = 0.106 to (±10,0) /xm. One beam is then kept 
fixed at this position while the other one is moved further 
outwards, dragging the associated vortex with it. At t = 500 
the former beam is ramped linearly down within At = 100 to 
release the single vortex close to the condensate center. The 
chosen grid size for this case is 749 x 749 resulting in a spatial 
discretization with spacing dx = 0.26. 


more or one fewer positively charged vortex compared 
with negatively charged vortices. To create such a non¬ 
neutral distribution in a controllable and repeatable man¬ 
ner we start by creating a neutral configuration as dis¬ 
cussed in the previous section and then take one vortex 
out of the condensate. We first demonstrate this proce¬ 
dure starting with a pair of oppositely charged vortices, 
removing one of the vortices, and leaving a single vortex 
that can then be re-positioned at will. The creation of 
a larger non-neutral configuration of vortices follows a 
similar principle. 

Figure [ 6 ] illustrates the ability to create a single vortex. 
First the protocol follows the generation of two vortices 
as in Fig. [2l When the beams have reached (±10, 0) /xm, 
one of the beams is kept stationary at its position while 
the other beam is moved further towards the edge of the 
condensate, dragging the trapped vortex with it. Ideally, 
the vortex would be dragged all the way out of the con¬ 
densate. However, shortly before reaching the edge of the 
condensate the vortex detaches from the beam and this 
edge vortex starts circulating indefinitely around the con¬ 
densate (see Fig. [ 6 ] where the vortex is barely visible in 
the last two density snapshots, but is clearly visible in the 
corresponding phase profile). This type of “detachment” 
is a general issue that is encountered when trying to “rid” 
the configuration of individual vortices by moving one of 
the beams out of the BEC in a direction that is normal to 
the edge of the BEC. A number of pointers about how to 
bypass this issue, most notably adjustments to the beam 
trajectory, and relying on the dissipation stemming from 
thermal excitations, is discussed below. As soon as the 
outgoing beam has reached the edge of the condensate 
the first beam is linearly and adiabatically ramped down 
to release the remaining single vortex to undergo free 
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FIG. 7: (Color online) Similar to Fig. [ 6 ] but in the presence 
of dissipation, as described by Eq. dS]) with 7 = 2 x 10“^. 


time evolution. 

While it is advantageous to reduce the beam veloc¬ 
ity when dragging the vortex out to allow the vortex 
to follow the beam trajectory closer to the edge of the 
condensate (at the cost of larger time scales), it is — 
typically, in our observations— not sufficient to drag the 
vortex all the way out of the BEC. This is due to the 
fact that the density becomes small at the condensate 
edge and thus the (density) contrast created by the laser 
beam close to this edge is too small to pin or drag a 
vortex. Nonetheless, we have found that the presence of 
dissipation (typically present in all BEC experiments), 
as described by Eq. (| 8 ]), can resolve this issue. Eigure[7| 
depicts the results corresponding to Eig. [ 6 ] with the ad¬ 
dition of dissipation with 7 = 2 x 10 “^. As it is evident 
from the figure, at t = 1200 the edge vortex is already 
very close to the edge of the condensate (as is visible in 
the phase profile) and has completely left the condensate 
by t = 1600. Eigure [HI shows a similar case but for the 
experimentally more realistic value 7 = 2 x 10 “^; see the 
relevant discussion in Ref. [33| , as well as in earlier works 
on coherent structures such as dark solitons in Ref. j^ . 
In this case, the vortex still is visible at t = 8000 and has 
completely left the condensate at t = 12000. We find the 
timescale for expelling the edge vortex from the conden¬ 
sate —with a starting point at about t = 600 where the 
dragging beam has reached the edge of the condensate— 
to be approximately a factor ten times larger than that 
for 7 = 2 X 10“^. As the physical picture is qualitatively 
the same but the (computational) timescales are consid¬ 
erable shorter in the following we will always opt to use 
the value 7 = 2 x 10 “^ to obtain results in the presence 
of dissipation. 

A more efficient protocol to drag a vortex out of the 
BEC relies on adjusting the laser beam path such that 
its trajectory becomes more azimuthal as the beam gets 
closer to the edge of the BEC and thus mimicking the 
natural tendency of the vortex to precess about the trap 
center. Eigure [9] depicts this improved protocol where 
one of the beams is moved on a circular path of diameter 
equal to the Thomas-Eermi radius (see white line) such 
that it leaves the condensate tangentially. Additionally, 
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FIG. 8: (Color online) Similar to Fig. [7| but in the presence 
of a more realistic dissipation value of 7 = 2 x 10“^ 



FIG. 9: (Color online) Improved protocol to create a single 
vortex. Same as Fig. [6] but with a circular trajectory (as 
indicated by the white line in the left panel) instead of a 
straight line to drag the vortex out of the condensate. Addi¬ 
tionally, the velocity of the dragging beam is reduced by 50% 
after t = 110. Att = 1000, when the latter beam starts to 
leave the condensate, the other beam is ramped down linearly 
within At = 100. 


the velocity of the dragging beam is reduced by 50% after 
t = 110. This results in the successful creation of a single 
vortex at t = 1620 as the second vortex seems to be 
completely removed from the BEG (it is neither visible 
in the density nor in the phase plot). 

For the generation of three, five or seven vortices with 
a charge imbalance of one, the idea is the same: start 
with a neutral configuration with for four, six or eight 
vortices (see Fig.[5|) and displace one beam outward with 
the aim of removing a single vortex from the neutral con¬ 
figuration. More specifically, after an even number of vor¬ 
tices is created, all but one of the beams then are kept 
stationary at their respective locations while the remain¬ 
ing beam is moved further outward to ideally drag the 
trapped vortex all the way out of the condensate. This is 
followed by linearly ramping down the other beams and 
releasing the vortex configuration in order for it to per¬ 
form free (i.e., unaffected by the beams beyond this time) 
time evolution. The results for this proposed methodol¬ 
ogy are displayed in Fig. [TOl The case for initially four 
(and subsequently three) vortices is depicted in the top 


two of panels in Fig. [TOl This case is more successful 
than for the creation of a single vortex, i.e., here the re¬ 
maining vortex has nearly vanished in the background. 
The case corresponding to initially six (and subsequently 
five) vortices is depicted in the middle rows of panels 
in Fig. [TOl In this case, the edge vortex is still visible 
at t = 1200 but cannot be distinguished from the back¬ 
ground at t = 1600. Finally, the bottom rows of panels 
in Fig. [To] illustrate the results for initially eight (and 
subsequently seven) vortices. This case example is again 
fairly successful in the sense that a potential edge vortex 
cannot be distinguished from the background. 


C. Comparison with vortex imprinting method 

The dynamics of vortices can also be explored by im¬ 
printing vortex solutions onto the ground state BEG in 
the presence of only the harmonic confinement. Such a 
technique has been used successfully in generating co¬ 
herent structures [38|, although it should be noted that 
typically in these experiments only a phase imprinting is 
induced. The latter necessitates the “morphing” of the 
density around the imprinting spot into a vortex profile, 
a process which, in turn, generates considerable sound 
wave emission clearly visible in some of the case exam¬ 
ples of Refs. [38|. Here, when we refer to imprinting, 
we more accurately mean both phase and density engi¬ 
neering, or “implanting” a vortex in the system, namely 
imprinting the phase, concurrently with modulating the 
density in order to produce an “as nearly exact as possi¬ 
ble” vortex waveform in the system. While the latter sce¬ 
nario is quite idealized and not straightforwardly achiev¬ 
able practically in the laboratory, our aim is to compare 
the dynamics of our produced vortices via the chopsticks 
method (containing the density modulations induced by 
the light beams etc.) to “target” vortex dynamics for a 
similar set of initial vortex locations. 

In the following we compare the time evolution that is 
initialized by the chopsticks method and compare it with 
the results obtained for the vortex implanting method. 
It should be stressed that while the chopsticks method 
generates (sometimes large) sound waves in the conden¬ 
sate as well as allowing vortices to move prior to fully 
ramping off the laser beams, the implanting method gen¬ 
erates the most pristine setting with minimal sound cre¬ 
ation and nearly “pure” vortex dynamics. Note that this 
is a numerical comparison only, as there have not been 
any demonstrations of phase engineering and imprinting 
of arbitrary vortex distributions into a BEG. We choose 
the case of three vortices, as displayed in the top rows of 
Eig. [TTJ To initialize the vortex implanting method we 
determine the positions and charge of the three vortices 
(resulting from two sets of split beams as explained in 
the previous subsection) at t = 1400. The ground state 
solution of the BEG is then multiplied with the corre¬ 
sponding (normalized) vortex solutions (obtained for a 
homogeneous background BEG) at the desired locations 





8 


nmnntTn 

rgm 






FIG. 11: (Color online) Comparison of the chopsticks proto¬ 
col and the vortex implanting method for the example of the 
three vortices from Fig. [T0]for the configuration at t = 1400 
after releasing the vortices from the beams (first column). For 
the vortex implanting method the time evolution is initial¬ 
ized by a wave function with two negatively charged vortices 
of charge located at (-10.025,-20.075) and (9.25,21.5) and 
one positively charged vortex located at (-8.45,21.5) (second 
column). Considerable differences in the time-evolution can 
be observed in comparison to the time-evolution generated by 
the chopsticks method and the vortex implanting method’s 
dynamics (columns three and four, respectively) due to the 
generation of sound waves in the former method. 
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FIG. 10: (Color online) Controlled generation of an odd num¬ 
ber of three to seven vortices for the same parameters and ini¬ 
tial configurations as in Fig. [5] Upper rows: At t = 80, three 
beams are kept fixed, while beam number four is moved fur¬ 
ther outwards dragging one vortex out of the condensate. The 
beam velocity is reduced to 0.25% of its original value. After 
linearly ramping down from t = 1200 to t = 1400 the other 
vortices are released from their beams. Middle rows: Simi¬ 
lar but for six vortices stemming from three pairs of beams. 
Five beams are kept fixed at t = 65 and linearly ramped down 
from t = 600 to t = 800, leading to a five vortex configuration, 
while the other is moved out of the condensate, again with 
reduced speed of 0.25% of the original value. Bottom rows: 
Similar to the previous cases but for four pairs of beams. 
Seven beams are kept fixed at t = 75 and linearly ramped 
down from t = 800 to t = 1000 while one vortex is dragged 
all the way out of the condensate, endowing the system with 
a configuration of seven vortices of alternating charge. 

(this is the “implanting” part). Fig. [TT] clearly shows that 
the two methods generate a non-trivially differing time 
evolution. While the qualitative structure of the vortex 
positioning and trajectories when using the chopsticks 
has been found to follow that of the “implanted” vortex 
configuration with the same initial vortex positions, the 


FIG. 12: (Color online) Vortex trajectories corresponding to 
the cases depicted in (a) Fig. [TT] and (b) Fig.[T3l The orbits in 
blue correspond to the chopsticks method while the red orbits 
correspond to the implanting method. 


density modulations induced by the vortex generating 
beams definitely affect the precise aspects of the vortex 
dynamics. The comparison of the corresponding trajec¬ 
tories using both methods is depicted in Fig. [T 2 ja). 

Next, and for reasons of completeness, we investigate 
the influence of dissipation on both cases. While the 
presence of dissipation with 7 = 2 x 10 “^ again helps to 
expel the edge vortex from the condensate that is visible 
at t = 1400 (and barely visible at t = 2000) for the 
vortex dragging method it does not have a big influence 
on the motion of the other three vortices closer to the 
BEC center. Also for the vortex implanting case the 
influence is negligible, cf. Fig. [131 Hence, unfortunately, 
in this case the presence of thermal excitations does not 
considerably alleviate this discrepancy. The comparison 
of the corresponding trajectories using both methods is 
depicted in Fig. [T 2 fb). 
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FIG. 13: (Color online) Same as Fig. ITT] but in the presence 
of dissipation with 7 = 2 x 10 “^. 



FIG. 14: (Color online) Principle of repository beams: Two 
stationary large beams with bw — 15 /xm and beam height 2/x 
serve as repository beams to collect vortices of equal charge. 
Vortices are created in oppositely charged pairs by ramping 
up two chopstick beams with bw — ^ /xm at (0, —30) /xm (at 
t = 0 they are already present) and by moving the beams to 
the pinning sites at (±30, 30) /xm. When the beams reach 
the pinning sites they are ramped down while new chopstick 
beams are ramped up at (0, —30) /xm. Al t — 220 one vortex 
is present in each repository beam, at t = 1870 seven vortices 
are present in each pinning site. If a certain critical number of 
vortices per pinning site is exceeded vortices start leaking out 
of the repository beams, cf., panel for t = 1980 where a vortex 
has been expelled from the repository beam. For this value 
of the chemical potential (/x = 4), the grid size is increased to 
1199 X 1199, and dx — 0.23. 


D. Repositories 

Lastly, we explore the possibility of depositing several 
vortices in so-called repository beams. The latter can 
be important for various reasons in the form of persis¬ 
tent currents [HI, [HI , but also towards the monitoring of 
dynamics of large vortex clusters aggregating with a sin¬ 
gle sign of vorticity, so-called Onsager-Kraichnan conden¬ 
sates, also thermodynamically representing negative tem¬ 
perature states; see, e.g., the recent analysis of Ref. 0. 

As an example, we investigate a configuration where 
two stationary repository beams with large beam waist 
are located in the BEG, cf. Fig. [HI The idea is to cre¬ 


ate one pair of oppositely charged vortices after another 
and deposit the positively charged vortices in one of the 
repositories while the negatively charged vortices are de¬ 
posited in the other repository beam. In this way, persis¬ 
tent currents can be obtained [HI • The vortex generation 
portion of the sequence follows the principle illustrated in 
Fig. |2j When the chopsticks reach the repository beams 
their position is kept fixed and they are linearly ramped 
down while a new chopsticks pair is linearly ramped up 
and another pair of vortices is created and dragged to 
the repository beams. This procedure can be repeated 
until a critical number of vortices is trapped in each 
repository beam. As soon as a certain number of vor¬ 
tices —depending on the size of the repository beam— 
is exceeded, vortices start leaking out of the repositories. 
A natural constraint in that regard is that the reposi¬ 
tory beam width should be larger than the product of 
the number of vortices times their corresponding length 
scale, i.e., the healing length of the BEG. 

To create a large number of vortices within this pro¬ 
cedure with as little disturbance of the condensate as 
possible, it is advantageous to use a larger chemical po¬ 
tential. In the following, we choose /x = 4. The size of 
the chopsticks beams is again b^ = Q /xm. In Eig. [14] we 
display results for repository beams with 6^ = 15 /xm. 
The time evolution is initialized by using the stationary 
state with two chopsticks beams present at (0, —30) /xm 
and repository beams located at (±30, 30) /xm. Eor a 
duration At = 175 the chopstick beams are moved to¬ 
wards (±30, 30) /xm, creating and dragging the first pair 
of vortices with them. At the position of the repository 
beams the chopstick beams are linearly ramped down 
within At = 100 and the respective vortices are conse¬ 
quently deposited in the repository beam. While the old 
chopstick beams are ramped down a new chopstick pair 
is linearly ramped up at the same initial position as the 
former, again within At = 100. Then the whole proce¬ 
dure is repeated and the second vortex dipole pair is cre¬ 
ated, dragged along and deposited, respectively, within 
the repository beams. Eor the chosen size of the reposi¬ 
tory beams seven vortices can be placed into each beam. 
After this critical number has been reached, vortices start 
leaking out of the repository beams. At t = 1980 the first 
vortex has been expelled from the repository beam. 

To accommodate even more vortices within the repos¬ 
itories, we increase the chemical potential to /x = 5 and 
the size of the repository beams to 30 /xm. Eor these pa¬ 
rameters ten vortices can be deposited in each repository, 
cf. Eig. [151 Next, both the chopstick and the repository 
beams are linearly ramped down within At = 100 to re¬ 
lease the vortices to their free time evolution. Eigure[T6] 
shows the same protocol as in Eig. [15] but in the pres¬ 
ence of dissipation with 7 = 2x10“^. It can be seen 
that the introduction of dissipation is, in fact, beneficial: 
the significant disturbances (relevant sound waves) to the 
atomic cloud due to the density modulations induced by 
the repositories, that can be seen in Eig. [15] are consid¬ 
erably smoothed out in the dissipative case. 
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FIG. 15: (Color online) Similar as Fig. [T4| but with larger 
chemical potential /i = 5 and larger repository beams 
with bw = 30 /xm. In this way ten vortices can be trapped per 
repository beam. Within At = 233 the dragging beams are 
moved toward the repositories and ramped up/down within 
At = 100 . At t = 3230, the repository beams are linearly 
ramped down within At = 100 to release the vortices to 
the free time evolution. However, importantly a consider¬ 
able amount of sound emission is present in the dynamics due 
to the significant density modulations imposed by the repos¬ 
itory beams. For this value of the chemical potential (/x = 5) 
the grid size is further increased to 1499 x 1499, resulting 
in dx = 0 . 21 . 



FIG. 16: (Color online) Same as Fig.fTSlbut in the presence of 
dissipation with 7 = 2 x 10“^. Notice that the sound waves 
have been chiefly expelled from the BEG, enabling a more 
pristine monitoring of the vortex dynamics. 


This specific configuration with two repository beams 
can easily be modified to other interesting configurations. 
Note that a repository beam sitting on the edge of the 
condensate can serve as a location to trap edge vortices 
that can arise when a vortex is supposed to be taken out 
of the condensate as in Fig. [6] in a reproducible manner. 
However, note that our results demonstrate that the pres¬ 
ence of dissipation helps to get rid of the edge vortices 
while the vortices closer to the center are less affected by 
the latter. 


IV. CONCLUSIONS & FUTURE CHALLENGES 

In the present work, we have numerically explored a 
mechanism that enables the production and manipula¬ 
tion of multiple quantized vortices, essentially at will, 
inside an atomic Bose-Einstein condensate. The use of 
lasers as vortex “optical tweezers” in a judicious manner, 
so as to gradually create the phase profile associated with 
a pair of oppositely charged vortices, as well as to pin vor- 
ticity and create persistent current, enabled us to locate 
vortices at various positions within the BEC by moving 
the laser beams. We were subsequently able to transfer 
these vortices at will, including moving them outside the 
BEC or positioning them within a repository beam. Re¬ 
peating the process either with multiple optical beams 
or after the delivery of the first pair to the repository 
beams, starting the process anew, we were able to pro¬ 
duce arbitrary neutral and non-neutral vortex distribu¬ 
tions. Naturally, the process has a number of limitations, 
such as the emergence of density modulations due to the 
carrying beams (which also to some extent affect the vor¬ 
tex motions), or in some cases the difficulty of carrying 
individual vortices outside the condensate (due to the so- 
called “detachment” from the beam in regions of low den¬ 
sity near the condensate rims. With respect to most of 
these aspects, the contribution of thermal fluctuations in 
the context of the so-called dissipative Gross-Pitaevskii 
equation is beneficial, enabling the outward motion of 
vortices and also the partial reduction of sound waves. 
The method is also generally limited by the number of 
laser beams of finite width that can be located and moved 
within a BEC of finite radius. 

This technique creates a broad set of possibilities that 
are quite worthwhile for subsequent experimental and 
further numerical exploration. As far as we know, no 
other technique available in the literature is as versa¬ 
tile towards creating/engineering multiple vortices while 
selecting their charge and position distributions at will. 
Constructing and understanding the dynamics of such 
vortices and vortex clusters [dj, m, especially in the 
context of anisotropy where they can be robust, e.g., in 
collinear states [dj, Q exploring more systematically 
their stability is now a tractable and experimentally re¬ 
alizable topic. Further studies oriented towards devising 
methods of decreasing the generation or effects of residual 
sound waves so as to enable a more direct engineering of 
free-vortex states and initialization and study of vortex 
dynamics would certainly be desirable. Finally, extend¬ 
ing these types of methods to higher dimensional BECs 
in three dimensions in the case of vortex lines and vortex 
rings would be of particular interest in its own right. 
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